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Abstract 

3 . 

We introduce the concept of the "polarized" distance, which distinguishes the orthogonal 

states with different energies. We also give new inequalities for the known Hilbert-Schmidt 
• »— i . 

distance between neighbouring states and express this distance in terms of the quasiprobability 

H ' 

distributions and the normally ordered moments. Besides, we discuss the distance problem 
in the framework of the recently proposed "classical- like" formulation of quantum mechanics, 
based on the symplectic tomography scheme. The examples of the Fock, coherent, "Schrodinger 
cats," squeezed, phase, and thermal states are considered. 
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1 Introduction 



Last years, an increasing interest to the problem of distance between quantum states is observed. 
Different motivations of this activity can be found in such fields as quantum cryptography, quantum 
communications, or quantum computing. Here we discuss the topic mainly from the point of view 
of quantum optics. In view of recent impressive progress in creating and detecting various types of 
nonclassical states of light or cooled particles in electromagnetic traps, the problem of measures of 
distinguishability or closeness between different quantum states becomes actual. For example, in 
quantum optics, the Glauber coherent states [jjj 

|a)=exp(-M 2 /2)£^|n> (1.1) 



are considered frequently as reference states (|n) means the Fock state with the definite number of 
photons), so that the (pure) states different from ( |1.1| ) are called sometimes as nonclassical states. 
But what is the quantitative measure of the "nonclassicality?" The simplest option is to use the 
so-called Mandel's parameter, Q = n 2 /n — n — 1, which equals zero for all coherent states, since 
they have the Poissonian photon statistics. However, this parameter is adequate for a limited class 
of states. Consider, for instance, the even and odd coherent states introduced in Q 

\a;±) = (2 [l ±exp(-2|a| 2 )])" 1/2 (|a) ±\-a)). (1.2) 

In this case, Mandel's parameter equals = ±2|a| 2 / sinh(2|a| 2 ), and it shows distinctly the 

qualitative difference between the states \a), \a; +), and \a; —), but only for small values of \a\. If 
\a\ S> 1, then ~ 0, although the states \a; ±) are still quite different from the coherent state. 
Moreover, for generalized coherent states (||, [|] 



a) 



00 

exp (-|a| 2 /2) J2 —j=. ex P i i( P ( n )] l n ) ( 1 - 3 ) 
v ' „_ n Vn! 



n=0 

we have identically Q = for any function c/?(n), although the state |5) may be essentially different 
from the Glauber state For example, the choice (p{2k) = (mod27r), (p(2k + l) = —tt/2 (mod2-7r) 
gives the so-called Yurke-Stoler state || 

™/ 4 (\a)+i\ -a))/V2 (1.4) 



\oi)ys = e 

which is considered, equally with the even and odd states, as a representattive of a large family of 
"Schrodinger cat states." 
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The concept of distance gives a possibility to characterize more precisely the neighbourhood or 
similarity between the quantum states. However, the existing approaches (see section seem to 
suffer from certain drawbacks. Some of the available definitions of a distance are too complicated to 
perform concrete calculations. On the other hand, some consequences of the traditional approaches, 
being correct mathematically, contradict the physical intuition. For example, the known definitions 
yield the same, at once, maximum possible value of the distance between any two orthogonal pure 
states, whereas from the physical point of view, the distance between the first and the 100th 
Fock states seems to be much greater than that between, say, the 100th and the 101th states. The 
distance measures based on the density operators alone are not sensible to the difference in energies. 

In the present paper, we propose new measures which distinguish different orthogonal states 
and which are simple enough to perform the calculations, at least for the most important families 
of states used in quantum optics. In our approach, the distance depends not only on the density 
operators alone, but also on some extra fixed positively definite operator. Of course, following 
this way we meet the problem of the nonuniqueness in the choice of this additional "polarization" 
operator. Nonetheless, such a nonuniqueness seems not crucial in many physical applications, 
where the special role of some operators (like the Hamiltonian or the quantum number operator) 
is evident from the beginning. Another goal is to provide an analysis of the distance problem in 
terms of the quasiprobability distributions and in the framework of the "classical-like" formulation 
of quantum mechanics proposed recently in ||. 

The paper is organized as follows. In section ^ we give a review of the existing approaches 
to the quantum distance problem. In section |3[ we concentrate on the properties of the Hilbert- 
Schmidt distance (HSD) and we express it in terms of the quasiprobability functions and ordered 
moments. In section we propose several definitions of the energy-sensitive distance in terms of the 
statistical operators (density matrices). In section the distinctions between different definitions 
are illustrated by examples of the Fock, coherent, "Schrddinger cat," squeezed, phase, and thermal 
states. The "classical-like" distances between quantum states are considered in section ^. The last 
section contains brief conclusions. 
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2 Previous approaches to the quantum distance problem 



The distance between two objects a and b is denned usually as a scalar real function satisfying the 
following properties: 



(I) d(a, a) = 0, d(a, b) > 0, if a ^ b, 
(II) d(a,b) =d(b,a), 
(III) d(a, 6) + d(6, c) > d(a, c) . 



(2.1) 
(2.2) 
(2.3) 



The property (III) has a clear geometrical meaning as the triangle inequality, and it implies rather 
strong limitations on the possible choice of the function d(a, b). If the "objects" a and b are different 
pure quantum states, then the distance must be some functional written in terms of the Hilbert 
space vectors, \a) and \b), representing the states. One should remember, however, that the set 
of quantum states is in one-to-one correspondence not with the whole Hilbert space of the wave 
functions, but with its projective factor space, since the vectors \ip) and e^lip) describe the same 
state. All the requirements are satisfied, e.g., for the Fubiny- Study distance [@, ^, ||| 

o\ 1/2 

(2.4) 



(2.5) 



d (FS) (Vi,V> 2 ) = v^(i-KVil^)! 1 

(sometimes the factor is replaced by 1 or 2), although a slightly different definition 



<i {mil %i,V>2) = inf 



\ih)-e^\rh) = v^l-KVil^)!) 172 



is also possible [JHJ. Taking a one-parameter family of states ip(t) generated by the time evolution 
operator, one obtains, both from ( |2.4| ) and (2.5), the infinitesimal distance along the evolution 
curve in the projective Hilbert space 



ds 



^2 - 2\mmt + dm 2 « 2^1 - \m)m + <**)>!■ 



(2.6) 



The definition fl2,6|) was used in studies devoted to the geometrical aspects of the quantum evolution 
and generalizations of the time-energy uncertainty relations (9|, 1C, 11, 13, 14, 15, 16, o]. 
For a family of states ip(s) dependent on a continuous vector parameter s = (si, S2, • • • , s n ) G 
R n , one can introduce the Riemannian metrics according to \\ip(s + ds) — vp(s)\\ 2 = jijdsidsj and 
measure not the "shortest" distance ( |3.1[ ), but the distance along a geodesies on a curved manifold, 
which can be much greater than the "shortest" one. The concrete examples of the geometries on 
the manifolds corresponding to the most known continuous families of quantum states (namely, 
coherent, squeezed, and displaced states) were studied in detail in p|, pi pi, M, M, E3 . 



Wootters [23] proposed the distance between the pure states in the form of the angle between 
the corresponding rays in the Hilbert space S w ^ l^)) = cos -1 [(V^l^i)!- F° r infinitesimaly 

close states, the differential form of this distance coincides (up to a coefficient) with ( ]2.6f) [f24|]. 
Recently, the Wootters and Fubini-Study metrics were compared in p5| . 

Now let us turn to the mixed quantum states, described by positively definite statistical oper- 
ators p with the unit trace: Tr/5 = 1. The first definition of the distance between mixed states in 
the physical literature, perhaps, was given in [E(| 



^( JMG) (/5i; 



P2j 



sup 



Tr [pi-p 2 ]i 



(2.7) 



Restricting the family of the bounded operators A in this definition by the projection operators 
E = E 2 , one obtains an equivalent definition [27| 



d^ (p 1; 



P2j 



sup 

E 



Tr([p!-p 2 ]i 



2 llPi -P2II1 



(2. 



where 



— Yl Wi I j the summation being performed over all the eigenvalues A n of 
the operator A. Actually, the right-hand side of equation (^j|) was used by Hillery [^] as a starting 
point in his definition of the distance between a state p and a given family of "classical" states p c i 



as S = infp cJ \\p — pctWi- More sophisticated definitions of the distance were given, e.g., in [29 



However, they are so complicated from the point of view of calculations, that no explicit examples 
were considered. 

One of the most frequently cited in the physical literature definitions is the so-called Bures- 



Uhlmann distance (BU-distance) j3l], plq| . It has the form (see also j27], 33, 34, p5| ) 

(Puh) = (2 - 2Tr V / ^ 7 Wp) 1/2 , (2.9) 

where the operator defined as the positively semidefinite Hermitian operator satisfying 

the relation = P- This operator is unique. Although the right-hand side of ( |2.9| ) seems 

asymmetrical with respect to p\ and P2, actually d^ BU ^(pi,p2) = d^ BU ^(p2,pi) fl35|] . For pure 
quantum states p^ = \ifi)(tp\, the BU-distance coincides with the "minimal" distance (^) due to 
p5, = pip- If one of the states is pure, then 



the relations p^ 2 



S BU Hm(^p) = ^[i-4^\ 



1/2 



(2.10) 



However, the calculations are much more involved in the generic case of nondiagonal statistical 
operators, so that the explicit forms of the Bures-Uhlmann distance were found only for finite- 
dimensional NxN density matrices (especially, for N = 2 and N = 3) [ 
squeezed thermal states [EM p3] and displaced thermal states [40]. 



36, B7f and recently for 



3 Distances based on the Hilbert- Schmidt norm 



A simple expression for the distance between quantum states, enabling to perform calculations for 
the most important classes of states (at least in the problems of quantum optics), is based on the 



Hilbert-Schmidt norm ||^4|| 2 = yTr(^4t^4). The Hilbert-Schmidt distance (HSD) of two statistical 
operators p x and p 2 is defined as §, |, |§, H> 0, IH, @|| 



d {HS \pi,P2) 



pi — P2II2 = { Tr (Pi ~hf } 
Tr( / 5?)+Tr(^)-2Tr( / 5 1/ 5 2 ; 



1/2 



(3.1) 



In particular (we write simply d instead of d^ HS ^ in all cases when it does not lead to a confusion), 



d(|V)M,/5)= 1 + Tr (p 2 



1/2 



< \/2[l 



\P 



1/2 



(3.2) 



so the HSD fl3.1|) goes to the Fubini-Study distance (|2.4| ) in the special case of two pure states. The 
possible values of the Hilbert-Schmidt distance are restricted by < d(p x , p 2 ) < \/2, the maximum 
\/2 being reached for any pair of orthogonal (pure) states. 

In many cases, it is convenient to describe the quantum states with the aid of quasiprobability 
distributions, which can be written as special cases of the general Cahill- Glauber s-distribution 



W(a, s) = Tr pT(a, s) 



(3.3) 



where 



T(a, s) 



d 2 ( 



7T 



■ exp 



C* (d-a) + 1\C\< 



a, ( are complex numbers and a, a' are the boson annihilation and creation operators (in one 
dimension, for simplicity). The choice s = (with a = (q + ip)/V2) yields the Wigner function p!| 
W(q,p) = J (in exp (ipu) (q — u/2\p\q + u/2). For s = —1, we have the so-called Husimi-Kano or 
Q-function [^] W(a, —1) = Q{a) = (a\p\a), whereas in the case s = +1 we arrive at the Glauber- 
Sudarshan function P(a) = W(a, +1) which yields the "diagonal" representation of the statistical 
operator 47] p = J P(a)\a}(a\d 2 a/7r. Using (|3.3| ) one can write the Hilbert-Schmidt distance in 
terms of integrals over the phase space: 

dq dp 



2tt 
d 2 a 



[W 1 (q,p)-W 2 (q,p)Y 



[Qi(a) - Q 2 (a)} [Pi(a)-P 2 (a)] 

7T 

d 2 ad 2 /3 „i a 1 2 



7T 7T 



[P l {a)-P 2 {a)\ [I\{fi)-P 2 {fi)]. 



(3.4) 
(3.5) 
(3.6) 



If one knows (e.g., from experimental data) all normally ordered moments M^ k ' 1 ^ = Tr (d) k a l p^ . 



then the statistical operator p can be reconstructed as follows [Etf, 49, pO 



k=0 1=0 



a>k,i 



min{fe,Z} 

E 



(-iy\i-j)(k-j\ 



j=0 jW(k-m-jy- 

Using this formula one can write the Hilbert-Schmidt distance in the form of a series 

(-l) s+k+l sl 



(3.7) 



(3.8) 



where AM^''' = M^ k — M% ' ■ For example, in the case of the coherent state \a) one has 
M{k,i) = a *k a i an( } gp converges to the closed expression ([0|). 

An advantage of the Hilbert-Schmidt distance is that it permits to obtain simple inequalities 
for the distances between neighbouring states. Consider, for example, the distance between an 
arbitrary state p and the vacuum state |0)(0|. Using formula ( |3.2| ) and the identities J2{n\p\n) = 1, 
J2 n ( n \p\ n ) = n, one can write the following chain of relations: 

1/2 



d(p,|0)(0|) < [2(1 - (0\p\0))} 1/2 = 

1/2 



2 EH/31 



< 



n\p\n) 



n=l 



n=l 



'2n. 



(3.9) 



This inequality is useful if n <C 1. For an arbitrary reference Fock state |n)(n|, one can prove in a 
similar way the inequalities 



d(p,\n)(n\) < V2[(0\p\0) +n-n(n\p\n)] 1/2 



1/2 



(3.10) 
(3.11) 



d(p, \n)(n\) <\/2 a n + (n - n) 2 

where cr n = n 2 — (n) 2 is the variance of the number operator in the state p. 

In general, one can identify the quantum state not necessarily with the statistical operator p, 
but with any function of this operator f(p). As a consequence, a whole family of the modified 
Hilbert-Schmidt distances can be introduced according to the definition 

1/2 



A/(pi,A0 = ll/(/5i)-/(P2)|| 2 =(Tr{[/(p 1 )-/(p 2 )] 2 } 



(Tr [/ 2 (pi) 



+ Tr 



f(p 2 ) -2Tr[/(p 1 )/(p 2 )] 



1/2 



(3.12) 



For pure states, -distances coincide with the Fubini-Study distance ( |2.4|) for any reasonable 
function f(p). However, for mixed states the new distances are essentially different. For example, 



choosing f(p) = p l l 2 we obtain the distance 



/2 



1/2 



(3.13) 



which coincides with the Bures-Uhlmann distance ( |2.9| ) for any commuting operators pi and p 2 
(remember that the pure state projection operators IV'KV'I an d \<^>){^\ do not commute if ^ \ 
If one of the states is pure, then 

d(\iP)(iP\,p) = V2[l- W|V>)] 1/2 , (3.14) 

so the inequalities (j3^)-( |3~Tl"l) hold for the d-distance, as well. 

4 Energy-sensitive distance between quantum states 

The Hilbert-Schmidt distance between any states cannot exceed the limit value y/2. In principle, 
one could "stretch" the distance between remote states, introducing some monotonous function 
F(d) with the property F(y/2) = oo. But such a simple modification yields the same (although 
infinite) distance for any pair of orthogonal states. 

To distinguish orthogonal states with different sets of quantum numbers, we have to break the 
symmetry of the Hilbert space with respect to "rotations" of the basis, i.e., to fix some "direction" 
given by a positively definite Hermitian "reference" operator Z . However, we still want to use the 
advantage of the Hilbert-Schmidt norm. So, we define the "Z-polarized" distance as 

1/2 



dz {pi,P2) 

Another possible definition is 

dz{f>uh) = 



(P! 



Pi. 



Tr(z [p 



Tr[Z p\ + p% - pip 2 - pif>\ 



~ -i2 
! " P2\ 

1/2 



(4.1) 



Tr Z 



Tr (Z 



\ . . .1/2.1/2 
P1+P2- Pi P 2 



.1/2 
Pi 



.1/2 .1/2 
P2 Pi 



1/2' 



P2 
1/2 



1/2 



(4.2) 



Evidently, both the definitions satisfy all the axioms due to the properties of the Hilbert-Schmidt 
norm (since we simply apply this norm to the "scaled" operators Z l l 2 p or Z l l 2 p 1 / 2 ). In the special 
case of pure quantum states pi = \ipi){ipi\, we have 

4 av>i>, iv> 2 » = d\ (i^i), iv 2 » = + <v> 2 i%2> 



■l$i\Z\^)^i) - ^ 2 \Z\^i)(i>i\i, 2 ) . 



(4.3) 



If Z coincides with the unity operator, ( |4,3D goes to the Fubini-Study distance (gjj). 

A possibility of using some extra operators to define the distance was mentioned in study [^] 
whose authors considered the construction Tr (^A [p\ — f>2\j- However, it was rejected on the grounds 
of the unboundness, if all observables A are admitted (the authors of |2(| started from the rough 
definition: "Two states are close to each other if all the expectation values of observables are close 
to each other"). Here we fix the operator Z, depending on the concrete physical problem. 

In the case of quantum optics, a natural choice of Z is the quantum number operator 

A> = a) a. (4.4) 

Then the TV-distance between the Fock states |n) and \m) reads 

d N (\m), |n)) = (1 - 5 mn )y/m + n. (4.5) 

We see that dAr(|m), |0)) > djv(|n), |0)) if m > n, i.e., higher the energy, more is the distance 
from the ground state. Nonetheless, the iV-distance also does not seem to be ideal. Consider, for 
instance, two Fock states with m, n ^> 1. Then d^dm), \n)) ^> 1, even if \m — n\ ~ 1. Such a 
property of the distance ( [4.1[ ) does not agree completely with our intuition. This drawback can be 
removed, if we assume the following definition: 

D% (p x ,p 2 ) = Tr (ApZAp) - ^ (A ^2 ■ ( 46 ) 



where Ap = pi — pi. The right-hand side of Eq. fl4.6| ) is nonnegative, since it can be written as 

Dl = Tr(Apf((z^-(z^)) 2 ), (4.7) 



where the average value is defined as (Z) = Tr (^ApZAp^j /Tr (Ap) 2 . We shall cautiously name Dz 
as a quasidistance , since we have no proof of the triangle inequality for any states. Applying ( |4.6| ) 
with Z = N to the Fock states, we obtain 

D N (\n),\m)) = \y/n - y/m\ /V2. (4.8) 

This expression obviously satisfies the triangle inequality. Moreover, it is in agreement with the 
representation of the Fock states in the phase space as circles whose radii are proportional to the 
square root of the energy 52]. In such a case, the distance between the 100th and 101th states 



is less than that between the ground and the first excited states. 
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A disadvantage of the definition ([O]) is that it complicates significantly calculations for non- 
Fock states. In the case of coherent states, the calculations are simplified if one slightly modifies 
the definition of the quasidistance in the following way: 

bl (ft, ft) = Tr (A^aAp) - lTl ^f ■ (4-9) 

Then 

D a (|a), \0)) = -j=\a-P\y/l + exp (-\a - /3| 2 ) . (4.10) 

The right-hand side of Eq. ( 4.10] ) is a monotonous function of \a — (3\, increasing from \a — f3\ at 
\a — P\ <C 1 to \a — P\/V2 at \a — fl\ S> 1. Although we have no proof that the quasidistance D a 
satisfies the triangle inequality (|2.3| ) for all states, we can prove that the function ( 4,1C| ) satisfies 
this inequality for all values of a and (3. 



1/2 



5 Examples 

5.1 Coherent and Fock states 

For two coherent states \a) and \/3), one finds 

d(\a), |/3)) = V2 [l - exp (-|a - /3| 2 } 

If |a — /3| <C 1, then d(\a),\(3)) ~ \[2\a — (3\ is proportional to the geometric distance of the 
displacement parameters a and /? in the complex plane, but it goes to V2 when \a - (3\ » 1. The 
A^-distance (4.3) between the coherent states is given by 

1/2 



(5.1) 



d N (\a),\(3)) = [|a| 2 + |/?| 2 -2Re(/3*a)exp(-|a-/?| 2 ) 



(5.2) 



so djv(|a),|0}) > djv(]/3), |0)) if \a\ > |/3|. The A r -distance is equal to the geometrical distance 
\a — P\ in the complex plane of parameters, if Re (a/3*) = (i.e., for orthogonal directions in the 
complex plane). In Fig. 1, we plot the HS- and A^-distances between the Fock state \m) and the 
coherent state la) 



d^ HS \\a),\m)) =V2[1 



d^da), \m)) = m + \a\ 



2m 



1/2 



ml 



2 a 



2m 



(m- 1)! 



>l 2 



1/2 



(5.3) 



(5.4) 



as functions of the mean photon number |a| 2 for fixed values of m = 1, 2, 3. The HS-distance has a 
minimum at |a| 2 = m. For small values of |a| 2 , we have d( HS \\a), \m)) > d^ HS \\a), \n)) if m > n, 
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but this inequality changes its sign if \a\ 2 is sufficiently large. The N-distance has a minimum only 
for m = 1, and the m-dependence is monotonous for all values of |ct| 2 . 

5.2 Squeezed vacuum states 

The squeezed vacuum state |53| j depends on the complex parameter £ with |£| < 1 

1/4 ~ yffitf 



(5.5) 



n=0 



The HS-distance between the states and IC2) reads (see also @g) 

v / 2|Ci-C 2 | 



d(\(x),\<2)) 



(|i - CiCII [|i - CiCII + vCTPKR^P)' 



1/2' 



(5.6) 



For I <C 1 and IC2I <C 1, this is the geometric distance of the complex squeezing parameters. Using 
the parametrisation £ = tanhre i<ji , r > 0, we have a simplified formula in the case of <pi = 4>2- 

2sinh \ (n - r 2 ) 



rf(ICi),IC 2 )) 



A/cosh(ri - r 2 ) 

For r 2 = 0, (5.7) gives the distance between the vacuum state and the squeezed state 
The TV-distance can be expressed as 



(5.7) 



4(ICi),lC2» 



IC1I 2 , IC2I 2 



+ 



1-ICil 2 1-IC2I 2 

+ 2 ' CiC i! 2 ~ r R ;i,? a) - idi 2 ) (i - . 

1 1 — Ci C2 r 



(5.8) 



If |Ci, 2 | *^ 1> then (5.8) has the same limit as the "unpolarized" Hilbert-Schmidt distance (|5T 
d,N ~ d ~ |Ci — C2 1 • However, for large values of the squeezing parameter these two distances 
become completely different. For example, in the special case A<p = arg(j — arg£ 2 = we have 



instead of (5/7) the expression (tj = \(j 



d% (|Ci>, IC2)) = sinh 2 n + sinh 2 r 2 



2 sinh n sinh r 2 
cosh 2 (n — T 2 ) 



(5.9) 



and c^at (|C), |0)) = sinh r. 



5.3 "Schrodinger cat" states 



Now let us consider the family of the "Schrodinger cat" states 



\a;ip) = [2 1 + cos(/3exp(— 2\a\ ) J Ma) 



(5.10) 
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The special cases of this family are even states (ip = 0), odd states (92 = ir), and the Yurke-Stoler 



states (f = vr/2). A more general set of states \a; r, <p) ~ \a) + re l ^\ — a) was studied in | 55| . The 
square of the distance between the coherent and cat states with the same values of the parameter 
a equals 



/2/I \ 1 w 1 ~ exp(-4|a| , 

d 2 (a;<p), a)) = — ' ' , 2V 5.11 

1 + cos V9exp(— 2\a\ z ) 



For the distance from the vacuum state, we obtain 



2[l-exp(-|a| 2 )] 
1 + cos </?exp(— 2\a\ 



d 2 (|a; , |0)) = i : L L:^/^,^ 2V (5.12) 



whereas the distance between two states with the same parameter a but different values of phases 
if 1 and if% reads 

,2/ \ [1 - exp(-4[q[ 2 )] [1 - cos(y9i - if 2 )\ 

a {fl,f2) = jr- 7 , yrhz , , , , 2 x r 5.13 

[1 + cos pi exp(— 2|ap)j [1 + cos ip 2 exp(— 2|ap)J 

For \a\ 3> 1, we have d 2 (fi,f 2 ) ~ 2 sin 2 Qtpi — f 2 \ /2). 

The iV-distances between the same states have an extra factor \a\: 

,2 /1 ^ t n \\ 1 |2 1 - cos(^exp(-2|a| 2 ) 

dt ( a: if) , 0) = a — —, ^4tt, (5-14) 

N " ^ 1 " 11 1 + cos</?exp(-2|a| 2 )' v ; 

,2 ( \ _ |«| 2 [1 + exp(-4|a| 2 )] [1 - cos(y?i - f 2 )\ , - 

"at IVl > <P2) — u— j / i |2m Fi 1 7 ol i2\r (o.lbj 

[1 + cos 991 exp(— 2|ap)J [1 + cos if 2 exp(— 2|a| z )J 

Now we have djq m \/2\a\ sin (\<f% — p 2 \ /2) for \a\ S> 1. 

Equations ( |5.1l| )- (|5.15| ) clearly show that the YS-states are intermediate between even and odd 

ones. Moreover, we see that the distance between the YS and the odd states with the same \a\ is 

greater than that between the YS and the even states, and the YS-state is farther from the coherent 

state than the even state (whereas the Mandel parameter does not distinguish the coherent and YS 

states at all). This example demonstrates how the concept of distance helps to understand better 

the properties of different families of quantum states and the mutual relations between them. 

5.4 Coherent phase states 

As a further example, we consider the coherent phase states |56f l 

00 

| E ) = Vl -ee* ^e n \n), £L|e)=e|e), |e| < 1, (5.16) 



n=0 

where 



E- = ^ |n — l)(n\ = (aa^ a 



n=l 
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is the Susskind-Glogower phase operator [57] which can be considered to certain extent as a quantum 
analogue of the classical phase e ltf '. The HS distance between the states |ei) and |£ 2 ) is given by 

I — o 1 &2 1 

It is proportional to the geometric distance of the complex parameters e\ and £ 2 for |ei 2I ^ 1- F° r 
any |e| < 1, the distance from the vacuum state is simply d(\s), |0)) = \/2 |e|. At the same time, 
the d/v -distance is given by 

i g 12 \ £ 12 

c& (|ei>, |e 2 )) = - — + 



1-|£!| 2 1-\S 2 \ 2 

+ Jl-kil 2 ) (l-k 2 | 2 ) [| £l £ 2 | 2 -Re( gl ^)] 
[1- 2Re(eie^) + |£ie 2 | 2 ] 2 



In particular, djv(|e), |0)) = |e| (1 — |e 



5.5 Thermal states 



2N-1/2 



The pure quantum state ( 5.16| ) has the same probability distribution | (n|e) | 2 as the mixed thermal 



state described by the statistical operator 

p=- \n)(n\ (5.19) 

provided that one identifies the mean photon number n with |£| 2 / (1 — |e| 2 ) |5q| . Moreover, the 
state ( 5.16| ) arises naturally as an exact solution to some nonlinear modifications of the Schrodinger 



equation [fffil , so it can be named also a "pseudothermal state" [M]. Therefore it is interesting to 
compare the expressions ( 5.17| ) and ( 5.18[) for the distances between "pseudothermal" states with 



the analogous formulae for the true thermal states. 



The HS distance between two states (5.19) reads 



' y/{l + W) (1 + 2n 2 ) (1 + ni + n 2 ) 
Although it is proportional to the difference of the mean photon numbers, it goes to zero when 
ni j2 — > 00 and \n\ — n 2 | = const. The distance to the ground state equals 

(n, 0) = - n , (5.21) 

and it tends to 1 when n — > 00, i.e., to the value which is \/2 times less than the maximal possible 
Hilbert-Schmidt distance. These results become clear if one remembers that highly mixed states 
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are located, in a sense, deeply "inside" the Hilbert space, since the density operators form a convex 
set with the pure states contained in the boundary Nonetheless, being justified from the 

mathematical point of view, these properties do not agree completely with our physical intuition, 
because usually we think on highly mixed states as almost classical ones (all the coherence is lost), 
which must be far away from the intrinsically quantum vacuum state. In particular, it seems a 
little bit strange that high temperature states are closer to the ground state than any pure Fock 
state. 

Using the modified HS distance ( 3.13[ ) (which coincides with the Bures-Uhlmann distance in 
the case involved) we obtain 



1 



1 + m + n 2 



1/2 



In particular, the distance to the ground state equals 

S BU H%o) = 



'2n 



n 



1/2 



(5.22) 



(5.23) 



vT+l(l + vT" 

and it tends to the maximal possible value y/2 when n — > oo. It is interesting to compare this 
formula with the analogous one for the "pseudothermal" state ( |5.16| ), but written in terms of the 
mean photon number: 

^(W.o)-^. 

We see that the BU-distance for the mixed states is always a little bit less than the distance between 
the vacuum and the pure pseudothermal state with the same value of n, in agreement with the 
reasonings of the preceding paragraph. For ni j2 3> 1, ( 5,22j ) is simplified 

v2 IvW — v^il 



m,n 2 ) 



y/ni + n 2 

The square of the iV-distance between two thermal states Q5.19 ) reads 



d N (n 1 ,n 2 ) 



\ni -n 2 W{ l + ni+n 2 ) 2 + 2nin 2 (l + 2ni)(l + 2n 2 ) 



(5.24) 



(5.25) 



(1 + 2fEi)(l + 2n 2 )(l + ni + n 2 ) 
As well as for the HS distance, the high temperature states occur not very far from the ground 
state: 



djv(n,0) 



n 



1 



when n 



oo. 



l + 2n 2 

At the same time, using the modified A^-distance ( f4.2[ ) we obtain the expression 

, . 2 

i — — — — 

/ + ni)(l + n 2 ) + ^Jr^n 2 ~\ 



d N (ni,n 2 ) = m +n 2 - 2^n\ii 2 



1 + ni + n 2 



(5.26) 
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which yields djsrin, 0) = n 1//2 , as well as for pure states. Analyzing formula ( [5,18| ) for the iV-distance 
between the "pseudothermal" states, one can check that the right-hand side attains the minimum 
(for fixed absolute values | £"1,2 1 ) if Re ( £ * £ 2) = I ^1^2 1 - This minimal distance can be written in terms 
of 711,2 i n the form very similar to fl5.26| ), but the last factor has the exponent 3 instead of 2: 

3 

I • (5-27) 



d'N min ( 



|£2}) = ni + n 2 - 2y/nin 2 



ni)(l + n 2 ) + ^nin 2 



1 + ni + n 2 



Since the fraction inside the parentheses does not exceed 1 (this is a consequence of the inequality 
n\ + n 2 > 2^nin 2 ), we have rfw miT i(|gi ), \^)) > djv(^i, ^2) for any pair of pure and mixed states 
with the same mean photon numbers. Equations ( |5.26| ) and ( |5.27 ) can be simplified for 724.2 3> 1: 



d N {m,n 2 ) ~ ni + n 2 



8 (nin 2 ) 3/2 
(ni + n 2 f 



d Nm i n {\ £ i)A £ 2)) ^n 1 +n 2 



_ _ 16 (nin 2 ) 



(ni + n 2 ) c 



If also \n\ — n 2 \ <C ni,2 , then we obtain approximate expressions resembling formula (4.8) for the 
gwasidistance between the Fock states, but with different coefficients 

\/3|ni -n 2 \ 



d 2 N (ni,n 2 ) Fa a/3 



rfjVmin(| £ l}) l £ 2)) 



'n 2 



'Tlx 



>n 2 



2 |m — 72-2 I 



V n l + V n 2 

The dependences of different distances between the vacuum and thermal or "pseudothermal" states 
on the mean photon number n are shown in Fig. 2. The distances of the pure states are larger 
than analogous distances of the mixed states with the same mean photon numbers, excepting the 
case of the djy -distance, which is the same both for the thermal and the phase coherent states. We 
can conclude that the -distance seems to be the most adequate from the physical point of view 
(at least for thermal states). 



6 "Classical-like" quantum distances 

It is accepted that quantum states are described in terms of the wave functions (state vectors in the 
Hilbert space) or density matrices (statistical operators). However, these complex-valued objects 
have rather indirect relations to the results of measurements, which are expressed in terms of real 
positive probabilities. Recently, a new formulation of quantum mechanics in terms of positive 
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classical probability distributions was proposed || 61, [6^]. It is a natural consequence of the 
concepts of the so-called symplectic tomography developed in [63, 64|. 

Let us introduce the two-parameter family of quadrature operators X^ IU = fj,q + up, —oo< 
\i,v < oo, where q and p are the usual coordinate and momentum operators (in one dimension, for 
simplicity). It can be shown that the probability distribution w^ v (X) of the real eigenvalues of the 
Hermitian operator X^ v is given by the following integral transform of the Wigner function: 

«VP0 = / ^6( M + V p-X)W(q,p). (6.1) 

The reciprocal transform 

W(q,p) = — J dX dfx dv exp[i(X - fxq - vp)]w^ v {X) (6.2) 

enables to express any Wigner function (and, consequently, any density matrix) in terms of the 
positive marginal probability distributions w^ v (X) which can be obtained, in principle, directly 
from an experiment with the aid of the homodyne detection schemes. Consequently, the description 
in terms of the family of classical distributions w flu (X) is completely equivalent to the standard 
description in terms of the density matrix or the wave function. This fact is the basis of the 



"classical-like" formulation of quantum mechanics [g, 61, 32, 65, 66, p7[ |. In this formulation, every 
quantum state is described not by a single complex-valued function tp(x) or p(x,x'), but by an 
infinite set of classical positive probability distributions w^ u (X), — oo < /x, v < oo. For example, the 
Fock state of the harmonic oscillator \n) is described by the family of the marginal distributions 

w $ (X) = w$ (X) —Ml I - X ) , (6.3) 

where H n (z) is the Hermite polynomial, while the marginal distribution w${X) of the vacuum 
state reads 

tog) (X) = 1 exp (—^-s) ■ (6.4) 

Now, considering the quantum states described by two different sets of the marginal distributions 
Wpv (X) and w^l (X) we can define the "classical-like" distance between these states as 

V c ab = Jdpdv gi/i, u)dP ab (w$ , «$) , (6.5) 

where d c ah {v}0 ,w$^ is some classical distance between the distributions w$(X) and w$,(X). 
A positive weight function <?(/x, v) is introduced to ensure the convergence of the integral over 
\i y v. Evidently, if the "partial distance" d c ah (w$ } w^vj satisfies the triangle inequality for all fixed 
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values of fj,, u, this inequality remains valid after multiplying by the positive function g((j,, v) and 
the subsequent integration over fi, v. 



Let us consider, for example, the "Kakutani-Hellinger-Matusita distance" 69] between two 
real nonnegative distributions Pi (x) and P 2 (x) 

1/2 



dx JP, (x) 



P2{x) 



(6.6) 



Taking into account the normalization condition we arrive at the "classical-like" analogue of the 
Bures-Uhlmann distance 

1/2 



V% = V2 J dfidug(fi, v) 1 - J dX^Jw$(X)w$(X) 



(6.7) 



The "classical-like" analogue of the JMG-distance ( |2.8| ) is obtained if one chooses for d c ah the classical 
Kolmogorov distance |5S| 

&K.{P\M = J dx\P 1 (x)-P 2 (x)\. (6.8) 

To illustrate the new approach, let us consider the P w -distance (|T7|) between two coherent 
states I a) and \(3). Each of these states is described by means of the families of the marginal 



distributions like 



■. exp 



X - X a (fj,,u) 



fi 2 + v 2 



X a (fj,, v) = V^(//Rea + vlma) . 



(6.9) 



Introducing the polar coordinates in the /if plane, \i = Roostf, v = Rsiwd, we see that the 
V n -distance between the coherent states depends on \a — j3\ only: 



roc p2n f 

T%p = J RdRj d$g(R,0) |2- 2exp 



--|a-/3| 2 cos 2 (tf-</?) 



1/2 



(6.10) 



(here 99 is the phase of the complex number a — (3). It is convenient to choose the weight function 
g(R,$) independent on $ and to impose the condition g(R)RdR = 1. Then for close coherent 
states, we have = 4|a — 0\ if \a — /3| -C 1. When \a — 0\ — ► 00, the V n -distance tends to the 
constant value 2tt^/2. 

The integral over fj,, v can be calculated explicitly for classical-like distinguishability measures 
(DM) which are defined by the same formula fl6,5| ) but without imposing the requirement fl2.3| ) (the 
triangle inequality) on the function d^ b (wffl , w^l J . The distinguishability measures are frequently 
used in the classical statistics and the information theory fl68| . Their applications to quantum 
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mechanical problems were discussed recently in [ 70 , 71 1 . The most known examples of classical DM 



are the Bhattacharyya coefficient [68fl 



B(Pi,P 2 ) = -In J dx^ P l (x)P 2 {x) 



(6.11) 



and the Kullback-Liebler distinguishability measure [68] 

J{P X ,P 2 ) = Jdx [Pi (s) - P 2 (x)} In ^ 



P(x) 



(6.12) 



For coherent states, both these measures yield similar dependences on the parameters a and (3, 
which differ only in a scale factor (we assume the same weight function g(ji, v) as above): 



47t|q-/?| 2 . 



(6.13) 



These quantum DM are unbounded when \ot — (3\ — > oo, but they do not satisfy the triangle 
inequality. 



7 Conclusion 

Let us summarise the main results of the paper. We have obtained new inequalities for the Hilbert- 
Schmidt distance and its modifications, which can be used for evaluating the "degree of proximity" 
between close quantum states. We have given new expressions for the Hilbert-Schmidt distance 
in terms of quasiprobability distributions and in terms of the ordered moments. We have con- 
structed the distances which are sensitive to the energy of quantum states. These "iV-distances" 
are unlimited and they distinguish different orthogonal states. Besides, we have shown how the 
concept of distance can be introduced in the framework of the new "classicaHike" formulation of 
quantum mechanics in terms of positive probability distributions of the rotated (in the phase space) 
quadrature operators. 
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Figure 1: The dependences of the iV-distance (three upper curves) and the Hilbert-Schmidt distance 
(three lower curves) between the coherent state \a) and the Fock states \m) with m = 1, 2, 3, versus 
the mean photon number in the coherent state \a\ 2 . The order of curves from bottom to top (in 
the part of plot nearby the vertical axis): the lower curves correspond to m = 1 while the upper 
ones correspond to m = 3. 



Figure 2: Different distances between the vacuum and the thermal (mixed) and pseudothermal 
(pure phase coherent) states versus the mean photon number. The order of the curves in the 
right-hand side of the plot (from bottom to top): iV-distance for the thermal state; the Hilbert- 
Schmidt distance for the thermal state; the Bures-Uhlmann distance for the thermal state; the 
Hilbert-Schmidt distance for the pseudothermal state; iV-distance for the pseudothermal state (it 
coincides with the modified iV-distance djv for the thermal state in the case concerned). 
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